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Abstract 

Substantial collective flow is observed in collisions between large 
nuclei at high energy, as evidenced by single-particle transverse mo- 
mentum distributions and by azimuthal correlations among the pro- 
duced particles. The data are well-reproduced by perfect fluid dy- 
namics. In a separate development, calculation of the dimensionless 
ratio of shear viscosity 77 to entropy density s by Kovtun, Son and 
Starinets within AdS/CFT yields rj/s = l/47r, and they conjectured 
that this is a lower bound for any physical system. It is shown that 
the transition from hadrons to quarks and gluons has behavior similar 
to helium, nitrogen, and water at and near their phase transitions in 
the ratio rj/s. Therefore it is possible that experimental measurements 
can pinpoint the location of this transition or rapid crossover in QCD 
via the viscosity to entropy ratio in addition to and independently of 
the equation of state. 

A transition from a phase of hadrons to a phase of quarks and gluons with 
increasing temperature and/or baryon chemical potential has been studied 
theoretically for three decades [1]. Calculations with effective hadronic field 
theories and with perturbative QCD give consistent and intuitively reason- 
able estimates of where this transition occurs, but at the present time only 
lattice gauge theory calculations on finite lattices of finite size can yield quan- 
titatively accurate numbers. For two flavors of massless quarks the phase 
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transition is of second order. For three flavors of massless quarks the phase 
transition is of first order. For the real world with nonzero masses for up, 
down and strange quarks the answer is not definitively known yet, but the 
answer is likely to be a rapid crossover from one phase to another without a 
rigorous thermodynamic phase transition, at least for zero baryon chemical 
potential. Indeed, there may be a line of first order phase transition in the 
plane of temperature T versus baryon chemical potential /x starting from the 
chemical potential axis and terminating at some critical point in the T — /i 
plane [2]. The Relativistic Heavy Ion Collider (RHIC) at Brookhaven Na- 
tional Laboratory was constructed explicitly to create quark gluon plasma. 
After more than five years of operation, what have the experiments told us? 

One of the amazing discoveries of experimental measurements of gold on 
gold collisions at RHIC is the surprising amount of collective flow exhibited 
by the outgoing hadrons. Collective flow is evidenced in both the single- 
particle transverse momentum distribution [3] , commonly referred to as radial 
flow, and in the asymmetric azimuthal distribution around the beam axis [4], 
quantified by the functions vi{y,pT),V2{y,PT), ••• in the expansion 

where y is the rapidity and is the transverse momentum. The func- 
tion V2{y = 0,pt), in particular, was expected to be much smaller at RHIC 
than it was at the lower energies of the SPS (Super Proton Synchrotron) at 
CERN [5]; in fact, it is about twice as large. Various theoretical calculations 
[6] support the notion that collective flow is mostly generated early in the 
nucleus-nucleus collision and is present before partons coalesce or fragment 
into hadrons. Theoretical calculations including only two-body interactions 
between partons cannot generate sufficient flow to explain the observations 
unless partonic cross sections are artificially enhanced by more than an or- 
der of magnitude over perturbative QCD predictions [7]. Thus quark-gluon 
matter created in these collisions is strongly interacting, unlike the type of 
weakly interacting quark-gluon plasma expected to occur at very high tem- 
peratures on the basis of asymptotic freedom [8]. On the other hand, lattice 
QCD calculations yield an equation of state that differs from an ideal gas 
only by about 10% once the temperature exceeds 1.5Tc, where ^ 175 MeV 
is the critical or crossover temperature from quarks and gluons to hadrons 
[9]. Furthermore, perfect fluid dynamics with zero shear and bulk viscosities 
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reproduces the measurements of radial flow and f 2 very well up to transverse 
momenta of order 1.5 GeV/c [10]. Parametric fits to the transverse momen- 
tum spectra of hadrons, such as pions, kaons, and protons, result in average 
transverse fluid flow velocities of more than 50/ 

An amazing theoretical discovery was made by Kovtun, Son and Starinets 
[13]. They showed that certain special field theories, special in the sense that 
they are dual to black branes in higher space-time dimensions, have the ratio 
r]/s — I/Att (in units with h — kB — c — 1) where 77 is the shear viscosity 
and s is the entropy density. The shear viscosity is rigorously given by the 
Kubo formula 



The connection between transport coefficients and gravity is intuitively clear 
since both involve (commutators of) the stress-energy-momentum tensor T*-^. 
They conjectured that all substances have this value as a lower limit, and 
gave as examples helium, nitrogen, and water at pressures of 0.1 MPa, 10 
MPa, and 100 MPa, respectively. Is the RHIC data, represented by radial 
and elliptic flow, telling us that the created matter has a very small viscosity, 
that it is a perfect fluid! 

The relatively good agreement between perfect fluid calculations and ex- 
perimental data for hadrons of low to medium transverse momentum at RHIC 
suggests that the viscosity is small; however, it cannot be zero. Indeed, the 
calculations within AdS/CFT suggests that r) > s/(47r). It will now be ar- 
gued that sufficiently precise calculations and mcasTircments should allow for 
a determination of the ratio r^/s as a function of temperature, and that this 
ratio can pinpoint the location of the phase transition or rapid crossover from 
hadronic to quark and gluon matter. This is a different method than trying 
to infer the equation of state of QCD in the form of pressure P as a function 
of temperature T or energy density e. 

What happens in a gas as compared to a liquid? The simplest way to 
understand the general behavior was presented by Enskog, as summarized 
in [14]. Shear viscosity represents the ability to transport momentum. In 
the classical transport theory of gases r^/s ~ Tlf^eeV, where Zfree is the mean 
free path and v is the mean speed. For a dilute gas the mean free path is 
large, /free ~ l/na, with n the particle number density and a the cross section. 
Hence it is easy for a particle to carry momentum over great distances, leading 
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to a large viscosity. (This is the usual paradox, that a nearly ideal classical gas 
has a divergent viscosity.) In a liquid there are strong correlations between 
neighboring atoms or molecules. A liquid is homogeneous on a mesoscopic 
scale, but on a microscopic scale it is a mixture of clusters and voids. The 
action of pushing on one atom is translated to the next one and so on until 
a whole row of atoms moves to fill a void, thereby transporting momentum 
over a relatively large distance and producing a large viscosity. Reducing the 
temperature at fixed pressure reduces the density of voids, thereby increasing 
the viscosity. (This is commonly experienced with motor oil in the cold winter 
months.) The viscosity, normalized to the entropy, is expected to be the 
smallest at or near a critical temperature, corresponding to the most difficult 
condition to transport momentum. What do the atomic and molecular data 
show? 

In figures 1 and 2 are plotted the ratio rj/s versus temperature at three 
fixed pressures, one of them being the critical pressure (meaning that the 
curve passes through the critical point) and the other ones being larger and 
smaller, for helium and water. (Other substances, such as nitrogen, behave 
similarly.) The ratios were constructed with data obtained from the National 
Institute of Standards and Technology (NIST) [15]. (Care must be taken to 
absolutely normalize the entropy to zero at zero temperature; we did that 
using data from CODATA [16].) The important observation [17, 18] is that 
Tj/s has a minimum at the critical point where there is a cusp. At pressures 
below the critical pressure there is a discontinuity in 77/ s, and at pressures 
above it there is a broad smooth minimum. This is an empirical observation. 
Figure 3 shows a plot oi rj/s (in units of its minimum value at the critical 
point) versus T/Tc for helium, nitrogen, and water along their critical isobars. 
Unfortunately there is no obvious universal scaling law. 

Numerical many-body simulations have been performed for a variety of 
systems. The result of a calculation for a 2-dimensional Yukawa system in 
the liquid state is shown in figure 4 [19]. Here F = Q'^/aT is the Coulomb 
coupling parameter with Q the electric charge and a the Wigncr-Seitz radius. 
It displays explicitly the dominance of potential contributions at low tem- 
perature and the dominance of kinetic contributions at high temperature, in 
agreement with the intuition provided by Enskog. 

The energy-momentum tensor density for a perfect fluid (which does not 
imply that the matter is non-interacting) is T'^" = —Pg^^ + wu^^u^ . Here 
w = P + e = Tsis the local enthalpy density and is the local flow 
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velocity. Corrections to this expression are proportional to first derivatives 
of the local quantities whose coefficients are the shear viscosity r] and bulk 
viscosity (Thermal conductivity is neither relevant nor defined when all 
net conserved charges, such as electric charge and baryon number, are zero.) 
Explicit expressions are to be found in textbooks [12] and they are quite 
lengthy. There is an ambiguity to the meaning of flow velocity for relativistic 
dissipative systems. In the Eckart approach is the 4-velocity for baryon 
transport. In the Landau-Lifshitz approach it is the 4-velocity for energy 
transport. Neither approach is more correct than the other. Both yield the 
same divergence of the entropy density current. 



d^s'' = ^{diU^ + dju'-lSijdku'') +^{dku'') +^{dkT + Tukf (3) 



Here x is the coefficient of thermal conductivity. This expression embodies 
the non-decrease of entropy. 

Perfect fluid dynamics applies when the viscosities are small, or when the 
gradients are small, or both. The dispersion relations for the transverse and 
longitudinal (pressure) parts of the momentum density are 



dimension of length and v is the speed of sound. Since w — Ts, and since 
usually the bulk viscosity is small compared to the shear viscosity, the dimen- 
sionless ratio of (shear) viscosity to entropy (disorder) t^/s is a good way to 
characterize the intrinsic ability of a substance to relax towards equilibrium 
independent of the actual physical conditions (gradients of pressure, energy 
density, etc.). It is also a good way to compare very different substances. 

How does this relate to hadrons and quark-gluon plasma? In the low 
energy chiral limit for pions the cross section is proportional to s/f^, where 
s is the usual Mandelstam variable for invariant mass-squared and /^r is the 
pion decay constant. The thermally averaged cross section is (a) oc T"^ / f^, 
which leads to 77/s oc (/^/T)^. Explicit calculation gives [20] 



where = rj/w and Di 




^ _ 15 ft 



(5) 
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Thus the ratio t]/ s diverges as T — ^ 0. At the other extreme hes quark-gluon 
plasma. The parton cross section behaves as a oc g^/s. A first estimate yields 
■q/s (X Asymptotic freedom at one loop order gives g"^ oc l/ln(T/Ar) 

where A^- is proportional to the scale parameter Aqcd of QCD. Therefore r]/ s 
is an increasing function of T in the quark-gluon phase. As a consequence, rj/ s 
must have a minimum. Based on atomic and molecular data, this minimum 
should lie at the critical temperature if there is one, otherwise at or near the 
rapid crossover temperature. 

The most accurate and detailed calculation of the viscosity in the low 
temperature hadron phase was performed in [20]. The two-body interactions 
used went beyond the chiral approximation, and inchidcd intermediate reso- 
nances such as the p-meson. The results are displayed in figure 5, both two 
flavors (no kaons) and three flavors (with kaons). The qualitative behavior 
is the same as in eq. (5). The most accurate and detailed calculation of 
the viscosity in the high temperature quark-gluon phase was performed in 
[21]. They used perturbative QCD to calculate the full leading-order expres- 
sion, including summation of the Coulomb logarithms. For three flavors of 
massless quarks the result is 

!? = s-i^ (6) 

We used this together with the two-loop renormalization group expression 
for the running coupling 

with At = 30 MeV, which approximately corresponds to using an energy 
scale of 27rT and Ajg^ = 200 MeV. The result is also plotted in figure 5. 
These results imply a minimum in the neighborhood of the expected value 
of Tf. pa 190 MeV. Whether there is a discontinuity or a smooth crossover 
cannot be decided since both calculations are unreliable near Tc- 

It is interesting to ask what happens in the large Nc limit with g'^Nc 
held fixed [22]. In this limit, meson masses do not change very much but 
baryon masses scale proportional to Nc\ therefore, baryons may be neglected 
in comparison to mesons due to the Boltzmann factor. Since the meson 
spectrum is essentially unchanged with increasing Nc, so is the Hagedorn 
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temperature. The critical temperature to go from hadrons to quarks and 
gluons is very close to the Hagedorn temperature, so that Tc is not expected 
to change very much either. In the large N^, limit the meson-meson cross 
section scales as According to our earlier discussion on the classical 

theory of gases, this implies that the ratio ri/s in the hadronic phase scales 
as A^"^. This general result is obeyed by (5) since it is known that scales 
as Nc- The large limit of the viscosity in the quark and gluon phase may 
be inferred from the calculations of [21] to be 



\sJqgp igm,f\n{B/{gm,)) 

where A and B are known constants but dependent upon Nf/Nc- Thus the 
ratio rj/s has a finite large Nc limit in the quark and gluon phase. Therefore, 
wc conclude that rj/s has a discontinuity proportional to if Nc — > oo. 
This jump is in the opposite direction to that in figure 5. 

So far the only quantitative results for viscosity in lattice gauge theory 
have been reported by Nakamura and Sakai [23] for pure SU(3) without 
quarks. This bold effort obtained r]/s ^ 1/2 in the temperature range 1.6 < 
T/Tc < 2.2, albeit with uncertainties of order 100%. Gelman, Shuryak and 
Zahcd [24] have modeled the dynamics of long wavelength modes of QCD at 
temperatures from Tc to 1.5Tc as a classical, nonrelativistic gas of massive 
quasi-particles with color charges. They obtained a ratio oi rj/s ~ 0.34 in 
this temperature range. 

Another interesting approach has been the computation oi rj/s in J\f — 
4 supersymmctric SU{Nc) Yang-Mills theory (SYM). In that theory it is 
possible to do calculations in the large coupling limit [25] and in the weak 
coupling limit [26] , and without too much imagination it is possible to find a 
paramtrization that interpolates smoothly between the two limits. However, 
SYM has no renormalization group running coupling, no asymptotic freedom, 
and no thermodynamic phase transition. Since it has so many more degrees 
of freedom than QCD as possible scattering targets, its viscosity to entropy 
ratio is much smaller than QCD at high temperature when compared at the 
same value of the gauge coupling. However, when the theories are compared 
at the same value of the Debye screening mass they do agree reasonably well 
[26]. 

It ought to be possible to extract numerical values of the viscosity in 
heavy ion collisions via scaling violations to perfect fiuid fiow predictions. 
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The program is to solve relativistic viscous fluid equations, with appropri- 
ate initial conditions and with a hadron cascade afterburner, over a range 
of beam energies and nuclei and extract r){T)/s{T) from comparison with 
data. This program is analogous to what was accomplished at lower energies 
of 30 to 1000 MeV per nucleon beam energies in the lab frame. At those 
energies, it was possible to infer the compressibility of nuclear matter and 
the momentum-dependence of the nuclear optical potential via the trans- 
verse momentum distribution relative to the reaction plane [28] and via the 
balance between attractive and repulsive scattering [29]. 

At RHIC some of the spcciflc proposals to extract or infer the viscosity 
to entropy ratio from data include: elliptic flow, Hanbury-Brown and Twiss 
(HBT) interferometry, single particle momentum spectra, and momentum 
fluctuations. Other possibilities which have not been worked on yet include 
jet quenching and photon and dilepton spectra. Some of the complications 
include the possibility that gradients are so large that the second-order dis- 
sipative equations of Israel and Stewart are necessary and that turbulence in 
the plasma may lead to an anomalous viscosity. Clearly this is an interesting 
and challenging goal but worth the effort! 

This work was supported by the US Department of Energy under grant 
DE-FG02-87ER40328. 
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Figure 1: The ratio 77/s as a function of T for helium with s normahzed such 
that s{T = 0) = 0. The curves correspond to fixed pressures, one of them 
being the critical pressure, and the others being greater (1 MPa) and the 
other smaller (0.1 MPa). Below the critical pressure there is a jump in the 
ratio, and above the critical pressure there is only a broad minimum. They 
were constructed using data from NIST and CODATA. 
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Figure 2: The ratio rj/s as a function of T for water with s normahzed such 
that s{T = 0) = 0. The curves correspond to fixed pressures, one of them 
being the critical pressure, and the others being greater (100 MPa) and the 
other smaller (10 MPa). Below the critical pressure there is a jump in the 
ratio, and above the critical pressure there is only a broad minimum. They 
were constructed using data from NIST and CODATA. 
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Figure 3: The ratio r\l s'm. units of its value at the critical point as a function 
of T jTc for three common systems. The curves are isobars at the critical 
pressure. They were constructed using data from NIST and CODATA. 
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Figure 4: The shear viscosity versus Coulomb parameter for a two dimen- 
sional Yukawa liquid (dusty plasma) from simulations by Liu and Goree [19]. 
Potential energy dominates at low temperatures and kinetic energy domi- 
nates at high temperatures. 
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Figure 5: The ratio r}/ s for the low temperature hadronic phase and for the 
high temperature quark-gluon phase. Neither calculation is very reliable in 
the vicinity of the critical or rapid crossover temperature. 
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